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ABSTRACT 

Young stellar systems are known to undergo outbursts, where the star experiences an increased 
accretion rate, and the system's luminosity increases accordingly. The archetype is the FU 
Orionis (FU Ori) outburst, where the accretion rate can increase by three orders of magnitude 
(and the brightness of the system by five magnitudes). The cause appears to be instability in 
the circumstellar disc, but there is currently some debate as to the nature of this instability 
(e.g. thermal, gravitational, magneto-rotational). 

This paper details high resolution Smoothed Particle Hydrodynamics (SPH) simulations 
that were carried out to investigate the influence of stellar encounters on disc dynamics. Star- 
star encounters (where the primary has a self-gravitating, marginally stable protostellar disc) 
were simulated with various orbital parameters to investigate the resulting disc structure and 
dynamics. Crucially, the simulations include the effects of radiative transfer to realistically 
model the resulting thermodynamics. 

Our results show that the accretion history and luminosity of the system during the en- 
counter displays many of the features of outburst phenomena. In particular, the magnitudes 
and decay times seen are comparable to those of FU Ori. There are two caveats to this asser- 
tion: the first is that these events are not expected to occur frequently enough to explain all FU 
Ori or EX Lupi; the second is that the inner discs of these simulations are subject to numerical 
viscosity, which will act to reduce the accretion rate (although it has less of an effect on the 
total mass accreted). In short, these results cannot rule out binary interactions as a potential 
source of some FU Ori-esque outbursts. 
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1 INTRODUCTION 

The traditional picture of star formation des cribes the free-f all col- 
lapse of a protostellar molecular cloud (e.g. lShu et al Jl 1987h . Con- 
servation of angular momentum will in general ensure that, within 
typical free-fall times of ~ 10 5 yr, the collapse produces a pro- 
tostar with protostellar disc. These formation rates are consistent 
with the observ ed statistics of prot ostellar objects, for example 
those in Taurus dKenvon et al.ll 199dh and with numerica l simula- 
tions ( iBate et al]2003l ; IStamatellos et al.l2007l : lBatel2009h . If these 
formation timescales are converted to average mass accretion rates, 
then it appears that the standard picture will form a star at the rate 
of 10 _ Mq yr~ , whereas current observations suggests an av- 
erage infall rate of 1O~ 6 M0 yr" 1 or lower. This leads to the re- 
alisation that accretion rates in protos tars are not constant, which 
has been confirmed ob servationally ( Herbigl ["l977l ; lArmitage et all 
1200 ll : IZhu et af]|2009i ). and that short periods of increased accre- 
tion, (accompanied by periods of mass pile-up where the infalling 
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matter is not accreted) can solve the apparent inconsistencies be- 
tween observation and theory. 

This also provides an explanation for the FU Orionis (FU 
Ori) outburst objects, Class to Class II protostellar objects which 
undergo a characteristic rapid rise in luminosity spanning up to 
five m agnitudes, and then decay on timescales of a few hundred 
years (Hartmann & Kenvon ll99rj) . Calculated accr etion rates for 
FU Ori objects are typically 1O _4 M0 j/r _1 or more dHerbid 19771 : 
iHartmann & Kenvonll996h, which show a strong increa se over typ- 
ical infall rates dKenvon et alJll993l : lFurlan et alj2008h . 

The precise origin of these FU Ori outbursts is not known, 
although almost all theories involve a protostellar disc, as 
these discs are expected to be present around the majority 
of e arly-type stars. The theorie s rang e from thermal instabil- 
ity dLinet ail [T9 85: Bell & Linl 1 1994), gravitational instabi lity 
dVorobvov & Basull200l |200d 120081 : iBolev & Durisenl l2008h . a 
combinat ion of gravitational in s tability and mag netorotational in- 
stability dArmitage eFaT] 1200 ll : Izhu et al.l [2009). even the pres- 
ence of a planet controlling accretion o nto the central star 
dClarke & Sverlll99d iLodato & Clarke! 1 2004h . All these theories 
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agree that triggering a disc instability event is crucial to provid- 
ing the observed accretion rates that cause the outburst. The aim of 
this paper is not to identify the correct theory of FU Ori objects, 
but to comment on what was once considered a potential cause of 
FU Ori: the encoun ter of a primary star plus protostellar disc with a 
discless secondary dKenvon et al . 1988; Bonnell & Bastien 1992b . 

Observations are beginning to reveal that outburst phenom- 
ena appear to belong to differe nt classes, e.g. FU Ori (FUors), 
EX Lupi (EXors) (Herbig 2007), and others. This work proposes 
that stellar encounters could produce outbursts that share many ob- 
servational feat ures with FUors and EXors (as was suggested by 
IPfalzned J2008h . combining treecode simulations with cluster dy- 
namics), and yet may belong to a different class. Results are pre- 
sented from a series of high resolution Smoothed Par ticle Hydro- 
dynam ics (SPH) simulations with radiative transfer dForgan et al.l 
2009) of the encounter of a star-disc system with a discless com- 
panion. The luminosity of these systems during the encounter is 
studied, and their potential for observation as a separate subclass of 
eruptive variable is discussed. 



2 METHOD 

2.1 SPH and the Hybrid Radiative Transfer Approximation 

Smoothed Particle Hydr o dynamics (SPH ) dLucvl 1 19771 ; 
iGingold & Monaghanl 1 1971 iMonagh an 1992) is a Lagrangian 
formalism that represents a fluid by a distribution of particles. Each 
particle is assigned a mass, position, internal energy and velocity: 
state variables such as density a nd pressure can then be c alculated 
by interpolation - see reviews bv lMonaghad d 19921. [2005). In these 
simulations, the gas is modelled by 500,000 SPH particles: the 
primary star (and the secondary companion) are represented by 
point mass particles, which ca n accrete gas par ticles if they are 
sufficiently close and are bound ( iBateetal .11 19951) . 

The S PH code used in this work is based on the SPH code de- 
veloped by |Bateetai] < ll995l) . It uses individual particle timesteps, 
and individually variable smoothing lengths hi such that the num- 
ber of nearest neighbours for each particle is 50±20. T he code uses 
a hyb rid method of approximate radiative transfer (For gan et al.l 
2009), which is built on two pre-existing rad iative algorithms: 
the po lytropic cooling approximation d evised by Stamatellos et al. 
d2007l), and flux-l imite d diffusion (e.g. IWhitehouse & Batell2004 
iMaver et al.ll2007l. see iForgan et"afl [2009 for details). This union 
allows the effects of both global cooling and radiative transport to 
be modelled without extra boundary conditions. 

The opacity and temperature of the gas is calculated using a 
non-trivial equation of state - this accounts for the effects of H2 dis- 
sociation, H° ionisation, He and He + ionisation, ice evaporation, 
dust sublimation, molecular absorp tion, bound-free and free-free 
transi t ions and electron scatte ring dBell & Linl ll994; Bole v et all 
120071 : IStamatellos et alj 12007b . Heating of the disc can also be 
achieved by P dV work, shocks and viscous dissipation. 

2.2 Initial Disc Conditions 

The discs used in this work were initially evolved in isolation for 
several Outer Rotation Periods (ORPs). This allows the disc to ap- 
proach an equilibrium state and becom e marginally unstabl e, and 
to develop steady-state spiral structures (Lodato & Rice 2004]). The 
discs extend from n n = 1 au to r ou t = 40 au. The disc used for 
most of the simulations described here has a mass of 0.1 Mq with 



a central primary star of mass 0.5 Mq. The initial surface density 
profile was chosen to be E oc r _1 , with a sound speed profile 
of c s oc r~3. One variant of the disc was also run with a disc 
mass of 0.2 A/q. These initial conditions (in particular the small 
disc radii) were motivated by recent work on disc fragmentation, 
which suggests that massive discs with radii of order ~ 100 au and 
greater will tend to fragment in t he outer regions, on timescales 
close to the dynamical t imesc ale dStamatellos et alj[2007l : IClarkd 
120091 ; iRice & Armitagd 1200^) . This is consistent with observa- 
tions that massive disc s tend to have outer radii less than 100 au 
dRodriguez et alj2005h . 



2.3 The Stimulus: Adding a Companion 

With the disc evolved into a quasi-steady state, a secondary star 
was then added (at a separation of at least 100 au to the primary, to 
prevent any non-linear perturbations in the disc by the secondary's 
sudden appearance). The secondary was added with several differ- 
ent initial conditions, such that the system exhibited different semi- 
major axes and eccentricities, comprising a suite of 4 simulations 
(see Table[T]for details). All encounters are prograde, and coplanar 
to the disc: all encounters are also parabolic (with the exception 
of simulation 4, which has an eccentricity of 7). This allows the 
characterisation of how outburst behaviour is dependent on initial 
orbital parameters. 



2.4 Resolution 

It is of utmost importance that SPH simulations involving mass ac- 
cretion appropriately resolve mass elements. To correctly model the 
accretion of matter onto both masses, the available feedstock for ac- 
cretion must be identified. The primary can accrete from the inner 
regions of the disc, which are the densest (and hence best repre- 
sented by SPH particles). In the innermost annulus (say up to 0.5 au 
from the disc's inner edge), there is 0.015 Mq of material (given a 
disc of mass 0. 1 Mq and a surface density profile E oc r _1 ). These 
simulations use 500,000 SPH particles: this mass therefore cor- 
responds to around 80,000 particles. The minimum mass element 
that SPH can resolve is typically one nearest neighbour group (e.g. 
lBate&Burkertlll997h . which in this work corresponds to 50 parti- 
cles. Therefore, the accretion feedstock for the primary constitutes 
around 1600 nearest neighbour groups, comfortably above the min- 
imum mass resolution. The problem of increasing numeri cal vis- 
cosity in the inner regions is unfortunately insurmountable dClarkd 
2009). The artifically hig h viscosity in the inner regions will pre- 
vent mass piling up (e.g. IRice & Armi tage 2009), so it should be 
expected that the primary's accretion rate will be underestimated 
(although the total mass accreted may not be affected, as the ac- 
cretion will begin earlier, and for a longer d uration). Comparing 
the co ntributions to the viscosity parameter a dShakura & Sunvaevl 
1973) from Reynolds stresses and from gravity shows that artificial 
viscosity dominates these discs within the inner 10 au. At R = 10 
au, CtReyn = a grav ~ 10~ 3 . 

The secondary can accrete from matter it encounters along its 
trajectory: the matter must come sufficiently close to become bound 
to the secondary before accretion is possible. This defines the sec- 
ondary's feedstock locale (approximately) as a semi-annulus in the 
disc, centred on periastron, with upper and lower radial bound- 
aries based on the secondary's gravitational influence. It is a semi- 
annulus because only disc material in the correct orbital phase will 
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Table 1. Summary of the orbital parameters investigated in this work. 



Simulation 
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Figure 1. Snapshots of the two discs used after several ORPs: the 0.1 Mq disc (top), and the 0.2 Mq disc (bottom). 



come into close proximity with the secondary. Material on the op- 
posite side of the disc during the encounter is typically teased into 
a tidal tail that is not accreted by either body. 

In more rigorous terms, a gas element i must satisfy the fol- 
lowing for capture by the secondary: 



„ „ , „ 1 2 GM 2 mi 

& = tikin, + tipot = 7rWlil> i2 

2 r i2 



<0, 



(1) 



where M 2 is the secondary mass, fi 2 and Vi 2 are the position and 
velocity of the gas element relative to the secondary. Assuming that 
the majority of capture occurs around periastron, then 



Vi 2 ~ \vn — v 2:Peri \ 

Ti 2 ~ |m — r 2 , peri\ 



(2) 
(3) 



where rn is the separation between the gas element and the pri- 
mary, and vn is the velocity of the gas relative to the primary. If 
the disc is Keplerian (and no radial motion is assumed), then 



vd = ruttn — 



GMi 



(4) 



The standard orbital equations can be used for v 2tPe ri- This de- 
fines the semi-annulus in the primary disc, where the secondary 



Duncan Forgan and Ken Rice 




Stellar Encounters in the Context of Outburst Phenomena 5 



exerts sufficient influence to potentially capture disc material. The 
upper and lower limits of this semi-annulus are found numerically 
by solving 



10" 



GMi 2GMi\ 



f2,peri I 



GM 2 

Ti2 



= 



(5) 



Having defined the semi-annulus, it is a simple matter to cal- 
culate its mass, given a surface density profile for the disc. In the 
case of a E oc r _1 disc with mass 0.1 Mq, this gives an avail- 
able mass of ~ 0.0057 Mq. Again, the simulations use 500,000 
SPH particles, and hence this corresponds to a particle number of ~ 
30,000. Therefore, the accretion feedstock for the secondary consti- 
tutes 600 nearest neighbour groups, again well above the minimum 
mass resolution. 



3 RESULTS 

3.1 Simulation 1 - The Reference Simulation 

Images of the reference simulation can be seen in Figure [2] As the 
seconda ry passes through the disc, it imparts sig nificant energy to 
the disc dLodato et al.l2007l ; lForgan & R ice 2009). It captures a sig- 
nificant secondary disc (roughly 0.006 Mq) during the encounter 
(slightly more than expected from the analysis given in the pre- 
vious section). This is due to the tidal forces exerted by the disc 
on the secondary as it reaches periastron, reducing its velocity and 
allowing more matter to be captured. The encounter draws out a 
significant tidal tail, and erases the spiral structure previously seen. 

The accretion rate of the primary and secondary can be seen 
in Figure[3] with corresponding accretion luminosities in Figure|4] 
The accretion luminosity is calculated using 



Lac 



1 GMM 

2 R 



(6) 



where R indicates the accretion radius of the object (taken for both 
objects to be 0.1 AU in these simulations). Particles are accreted if 
they pass within the accretion radius of the object, and are gravita- 
tionally bound. R is held constant throughout the simulation, only 
M and M change. As the luminosity varies linearly with the accre- 
tion rate, the accretion luminosities show similar peaks and troughs 
as the accretion rates (with the important modification that the pri- 
mary's larger mass makes it more luminous). 

As a guide, Figure [4] has several important events delineated 
by vertical lines: 

• t ~ 1209 yr - the secondary begins its ingress to the disc, es- 
tablishing a non-zero accretion rate. 

• t ~ 1256 yr - the secondary reaches periastron. The disc's tem- 
perature also peaks at this time. 

• t ~ 1279 yr - The secondary reaches its peak accretion rate: 
the matter that forms the secondary disc (located in the outer tidal 
tail generated at the location of the secondary) is in the process of 
infall onto its new parent star, not taking a well defined spheroidal 
shape until it has begun to exit the disc. The disc's luminosity has 
returned to near pre-encounter levels. 

• t ~ 1352 yr - the secondary begins its egress from the disc, 
and its gravitational influence diminishes. The primary disc must 
now begin readjustment - an increase in outward angular momen- 
tum transport gives a corresponding increase in inward mass flux, 
boosting the primary's accretion rate. The rate has a rise time of 
approximately 25 years. 




Figure 3. Accretion rates for the primary and secondary in Simulation 1, 
The black line denotes total mass accretion, the blue line indicates the pri- 
mary accretion, the red line indicates secondary accretion. 



• t > 1352 yr - the luminosities decay with timescales of hun- 
dreds of years. 

It must be emphasised that there are two separate accretion events, 
with distinct characteristics: the secondary's accretion undergoes 
enhancement for a period of around 150 years, reaching a maxi- 
mum accretion rate of 5 x Mqi/t^ 1 at periastron. Its accre- 
tion rate curve is a superposition of three characteristic features: 

(i) A smooth symmetric curve, width approximately 110 years. 
This is linked to the secondary's passage through the smooth com- 
ponent of the disc. 

(ii) A series of small spikes, widths of 1 year or less. These are 
due to the secondary's passage through overdensities caused by the 
disc's spiral structure. 

(iii) A major spike, width approximately 10 years. This is 
caused by the formation of the secondary disc. 

The primary's accretion curve consists of a steady increase 
to values of ~ 10~ 5 A/qj/t- -1 as the secondary leaves the disc, 
with a rise time of order ~ 10 years. This is followed by a slow 
decay with timescales in excess of 300 years. The accretion curves 
of the subsequent simulations discussed share some or all of these 
characteristic features. 

It should again be noted that these simulations are subject to 
numerical viscosity in their inner regions, preventing mass build up. 
Some theories suggest that this pile up is crucial to activating MRI , 
which causes the outburst dArmitage et alj|200ll : lzhu et alj|2009l) . 
The artificially high viscosity in the inner regions can be thought 
of as acting like MRI, which activates at the "wrong" temperatures, 
facilitating mass accretion without pile up, and hence underestimat- 
ing accretion rates during the outburst. These simulations can then 
be thought of as presenting the lower bounds of the accretion rate 
(in the limit where MRI is overactive). 

The tidal interactions between the disc and secondary cause the 
secondary to be captured on an eccentric orbit: this facilitates study 
of the repitition of this outburst event over several orbits. Figure|5] 
shows the evolution of the accretion rate over three periastra pas- 
sages; it can be seen that the magnitude of the outburst reduces sig- 
nificantly with each passage. The ability of the secondary to strip 
mass from the outer disc diminishes with each passage, partially 
due to the depletion of the outer disc itself, as well as the secondary 
disc's influence in regulating mass flow. Without mass stripping, 
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Figure 4. Accretion luminosities for the primary and secondary (and disc 
luminosity) in Simulation 1. The black line denotes total luminosity, the 
blue line indicates the primary accretion, the red line indicates secondary 
accretion, and the green line indicates the disc luminosity. 
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Figure 5. Long term evolution of the accretion rates for the primary and 
secondary in Simulation 1 . The black line denotes total mass accretion, the 
blue line indicates the primary accretion, the red line indicates secondary 
accretion. 



the secondary cannot accrete, and the disc is not obliged to signifi- 
cantly readjust its mass distribution, preventing the primary accre- 
tion event. This would suggest that these outbursts are not easily 
repeatable. 



3.2 Simulation 2 - Close Periastron 

Can repeatability be attained by reducing the secondary periastron? 
It is not unreasonable to assume that a closer approach allows the 
secondary to accrete from a more plentiful mass supply, potentially 
improving its ability to create repeatable outbursts. The secondary 
reaches periastron at t ~ 1 155 yr, resulting in the first spike in Fig- 
ure[6] Feature (ii) is less prominent, as the secondary passes through 
weaker spiral structure in the inner regions. Again, the peak accre- 
tion rate occurs when the secondary disc begins its infall (feature 
(iii)), giving a similar order of magnitude increase in accretion rate 
as Simulation 1. 

The second encounter is of similar magnitude, with an ac- 
cretion peak at t ~ 1691 yr. Note that the accretion rate peak is 
smoother than the first: the spiral structure in the primary disc has 




2000 
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Figure 6. Accretion rates for the primary and secondary in Simulation 2. 
The black line denotes total mass accretion, the blue line indicates the pri- 
mary accretion, the red line indicates secondary accretion. 



been almost completely erased, so the secondary sees a smooth dis- 
tribution of mass along its trajectory. Also, the secondary disc can 
maintain its shape during the encounter, reducing the infall onto 
the secondary. These facts combined eliminate the possibility of 
seeing features (ii) and (iii) in the second peak: all that remains 
is the smooth component (feature (i)). However, this magnitude of 
outburst cannot be maintained through subsequent orbits, decaying 
in similar fashion to Simulation 1 . Repeatability therefore seems to 
be limited for outbursts of this type. 



3.3 Simulation 3 - A More Massive Disc 

As the important factor is enhanced accretion, will adding more 
mass to the disc result in a larger accretion event? The results of 
Simulation 3 (where the disc is double the mass of Simulation 1) 
indicate some important differences (Figure[7]l. The disc's lack of 
high-m spiral modes compared to the lower mass disc (see Figure 
[]} prevents feature (ii) from being presented here. Periastron occurs 
at t ~ 1284 yr, where no significant accretion peak can be seen. The 
peak accretion rate is again seen when secondary disc infall oc- 
curs at t ~ 1338 yr. Overall, the features of this outburst event are 
similar in magnitude and duration to those seen in the simulations 
with a less massive disc, in particular the peak accretion rate of the 
secondary and its decay timescale. This would suggest that the be- 
haviour of the secondary during the outburst is relatively insensitive 
to disc mass (possibly because the secondary is already operating 
at peak accretion efficiency at lower disc masses). With that said, 
the pre-outburst and post-outburst accretion rates of the primary are 
slightly higher, and its decay timescale is slightly shorter. Depend- 
ing on the orientation of this system to the observer, this will have 
important implications for observation (see Discussion). 



3.4 Simulation 4 - A Hyperbolic Encounter 

The previous sections of this paper have indicated that accretion ef- 
ficiency of the secondary is linked to the velocity of the secondary 
relative to the disc. If the secondary moves through the disc too 
quickly, it may be expected that the accretion event will be reduced 
in magnitude in comparison to the other simulations shown. Simu- 
lation 4 was run to discover the effects of increased orbital velocity, 
by specifying a hyperbolic encounter (e = 7). Figure [8] indicates 
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Figure 7. Accretion rates for the primary and secondary in Simulation 3. 
The black line denotes total mass accretion, the blue line indicates the pri- 
mary accretion, the red line indicates secondary accretion. 
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Figure 8. Accretion rates for the primary and secondary in Simulation 4. 
The black line denotes total mass accretion, the blue line indicates the pri- 
mary accretion, the red line indicates secondary accretion. 



that the accretion efficiency is indeed reduced. The peak accretion 
rate is lower, and the duration of the accretion event is also reduced, 
lasting around 50 years. There is no evidence of feature (ii): the 
high velocity of the secondary's motion (and the limited time reso- 
lution of the data) prevent the detection of these peaks. There is no 
peak associated with the secondary's periastron (t ~ 1110 yr). 

The reduced peak is again associated with the infall of matter 
onto the secondary after the encounter (feature (iii)), but no disc is 
formed. The high-velocity encounter essentially destroys the disc, 
throwing significant amounts of material to large distances. It is 
this halo of matter which the secondary accretes from, but not ef- 
ficiently or for any length of time, as the velocity dispersion of the 
material is quite large. 

The primary's accretion behaviour remains similar to the other 
simulations, despite the dispersive action of the secondary. Indeed, 
the act of removing disc material may help in the detection of such 
an event (see Discussion). 



4 DISCUSSION 

4.1 The Potential for Observation 

4.1.1 Frequency of Occurence 

For outbursts from stellar encounters to be observed and cor- 
rectly classified, their occurence must be sufficiently frequent in 
an observer's field of view. To estimate the frequency of stel- 
lar encounters in a star cluster, the collision rate calculations of 
ICrarke&Pringld l ll99ll) are employed. They calculate the collision 
rate Vhu of a star-disc system with discless companions in a star 
cluster with stellar number density no, and a Gaussian velocity dis- 
tribution, characterised in one dimension by the velocity dispersion 



T hit = To 1 + 



GM, 



where 



r = 



4^Tvn GM t R d 



(7) 



(8) 



The discs used in these simulations has Rdisc = 40 au, and 
M* = 0.5Mq. To mimic the core of an open cluster, the free 
arameters are selected to be no = 100 pc~ 3 , V* — lfcms -1 



Binnev & Tremaind 1 19871 ; Idarke & Pringld Il99lh . This yields 
To = 3.01 x 1(T 4 Myr -1 , and Y hit = 3.28 x 10~ 4 Myr _1 . 
Assuming that the average protostellar disc has a lifetime of ~ 1 
Myr, Yhu is the probability that one disc will undergo an encounter 
in its lifetime. Therefore, out of roughly 3000 stars with discs, one 
will experience an encounter in its lifetime. This would imply that 
in an average open cluster, there will be at least one encounter- 
driven accretion event per Myr. Assuming that these encounters are 
randomly distributed in inclination, it should be expected that only 
around 6% of these encounters will be coplanaiQ. Finally, the dura- 
tion of these encounters must be taken into account: the primary's 
accretion rate remains enhanced for time intervals of At = 500 
years, which improves the chances of detection somewhat. If the 
probability of detecting an event is defined as 



Pob 



^events^t 1 — r hit stars finest 



(9) 



With fi nc = 0.06, and At = 500 yr, then for N stars = 3000, this 
yields P b s ~ 10 -4 . This shows that encounters of this type should 
not be frequently observed. This is of course an oversimplification: 
it does not account for a distribution of disc radii or masses, nor 
the subclustering that exists in bound systems. It also does not re- 
flect the young ages of most FU Ori systems (in general less than 
half the disc's lifetime). However, it provides a sufficient order-of- 
magnitude estimation to illustrate the rarity of these events. 

4.1.2 The Problems of Obscuration 

The accretion luminosities shown in the previous sections are in- 
trinsic luminosities: they do not account for the effects of optical 
depth. Consider the case where the disc is face-on to an observer: 
the primary star resides in a gap at the centre of the disc, and so 
the optical depth to the observer is relatively low. The secondary, 
however, penetrates the disc, accumulating matter and heating the 
surroundings. The Rosseland mean optical depth of the secondary 



For an encounter to be coplanar, the inclination must be lower than the 
disc's opening angle, defined by its aspect ratio ^ = This can be 
justified by rewriting equation |5]f or an inclined secondary 
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at periastron reaches r ~ 300. In this optically thick regime, the 
secondary is strongly obscured: any detectable emission will be re- 
processed and emitted by the disc at longer wavelengths. 

The effect of obscuration increases with inclination to the ob- 
server. The optical depth of the primary (and secondary) to the ob- 
server can increase to r ~ 10 at edge-on, completely shielding the 
accretion event from observers. This would suggest that encounters 
have to occur within a restricted range of inclinations to the ob- 
server in order to be directly observable. 



4.1.3 Indirect Observational Signatures 

If the secondary's accretion event is mostly screened by the disc, 
then what observational signatures can be identified? The disc SED 
for Simulation 1 was calculated by assuming that the disc emits a 
blackbody spectrum at each annulus, and integrating the contribu- 
tions: 

Fx = J 2nrB x (T eff (r))dr (10) 

The effective temperature of the spectrum is defined using 
midplane variables: 

Te„(r) = . r 7" W 1 . - (11) 

Tmid(r) + T mid (r) 

The effect of the encounter is to steepen the surface density 
profile, which in turn steepens the optical depth profile. This allows 
the cooler outer regions to radiate more efficiently, while the inner 
regions experience stronger screening. The result is a boost in the 
disc's flux at longer wavelengths, increasing the flux at 850/im (for 
example) by a factor of 10. Add to this the obscuration of the sec- 
ondary (especially its shorter wavelengths), and it appears that the 
outburst event is best searched for in the far-IR to sub-mm regimes. 

However, the effects of stellar irradiation are not included in 
the simulation. Can a substantial increase in reprocessed emission 
and scattered light be expected from the encounter? Stellar irradi- 
ation (and reprocessed emission) can be esti mated as a funct ion of 
temperature, assuming radiative equilibrium dlda & LirJl2004l) : 

T w- 280 (iiu)"(fe) 3 (12) 

In the initial stage, the primary is the sole contributor: the pri- 
mary's initial accretion luminosity (plus an intrinsic luminosity of 
a similar mass main-sequence star) gives 

T(r) « 300 (^f * (13) 

By comparison, the disc reaches temperatures of over 1000 K 
in the inner 3 AU initially (from dynamical heating sustained by 
marginal instability), so it is expected that reprocessed emission as 
a result of accretion should not be significant initially. During the 
encounter, the secondary's accretion luminosity increases to values 
of ~ 40L Q . This implies that stellar irradiation should be responsi- 
ble for temperatures of around 700 K at 1 AU from the secondary. 
Compressive heating increases the temperature in this region to a 
similar value: this suggests that reprocessed starlight will play a 
significant role, and should also be considered in terms of observa- 
tions. 



5 CONCLUSIONS 

This paper has investigated the possibility of a stellar encounter 
(where one participant has a protostellar disc) being the progenitor 
of an outburst phenomenon. The outbursts described here have sev- 
eral key features in common, originating from two distinct accre- 
tion events (for the primary and secondary respectively), indepen- 
dent of the secondary's orbital parameters. The secondary's accre- 
tion rate grows and fluctuates as it traverses the disc's spiral struc- 
ture: the peak occurs when the secondary disc is being formed from 
infalling stripped primary disc material, and the accretion rate then 
decays over several hundred years. The primary accretion rate (per- 
haps the easiest to observe) increases slowly with a rise timescale 
of tens of years as the primary disc readjusts to mass stripping, and 
decays on a longer timescale than the secondary. 

It has been established that these stellar encounters can en- 
hance accretion rates to levels corresponding to FU Ori and EX 
Lupi outbursts, and that they have observational signatures, that al- 
though low in probability to detect, are nonetheless possible in prin- 
ciple to see. But, it is important to emphasise that these encounters 
cannot be responsible for all FU Ori phenomena (or EX Lupi for 
that matter), for the following reasons: 

(i) This type of stellar encounter is too infrequent to explain the 
catalogue of outbursts currently observed, 

(ii) Such encounters cannot maintain the rapid periodicity or re- 
peatability required of some outbursts without significant decay of 
the outburst strength, 

(iii) This origin would predict the detection of a companion 
(perhaps in the infrared or submillimetre) for every outburst host. 
This is not the case; out of at least twe nty FUors, only seven have 
a confirmed companion (Pfalzner 2008). 

Despite this, outbursts from stellar encounters can mimic FU Ori 
well, with the correct general behavioural trends. They may also 
potentially have the same trigg ering mechanism - mass pile u p 
leading to MRI activation ( Armitage et al. 2001; Zhu et al J l2009h . 
although higher resolution simulations (which resolve the inner re- 
gion and reduce the artificial viscosity) are required to confirm this. 
Taking the fraction of FUors with a companion as a guide, it is 
estimated that at most 30% of outbursts detected could be due to 
an encounter or binary (with the actual figure presumably much 
smaller). The outbursts identified here represent a subtly different 
type of object, that although not currently detected, may be detected 
in future large-scale surveys of star-forming regions at far-infrared 
and sub-mm wavelength. 



ACKNOWLEDGEMENTS 

All simulations presented in this work were carried out using 
high performance computing funded by the Scottish Universities 
Physics Alliance (SU PA). Surface density plots were made using 
SPLASH jPriceil2007h . The authors would like to thank Philip Ar- 
mitage for useful discussions which helped to refine this work. 



REFERENCES 

Armitage P. J., Livio M., Pringle J. E., 2001, MNRAS, 324, 705 
Bate M. R., 2009, MNRAS, 392, 1363 

Bate M. R., Bonnell I. A., Bromm V., 2003, MNRAS, 339, 577 
Bate M. R., Bonnell I. A., Price N. M., 1995, MNRAS, 277, 362 
Bate M. R., Burkert A., 1997, MNRAS, 288, 1060 



Stellar Encounters in the Context of Outburst Phenomena 9 



Bell K. R., Lin D. N. C, 1994, ApJ, 427, 987 

Binney J., Tremaine S., 1987, Galactic Dynamics. Princeton Uni- 
versity Press 

Boley A. C, Durisen R. H., 2008, ApJ, 685, 1193 

Boley A. C, Hartquist T. W., Durisen R. H., Michael S., 2007, 
ApJL, 656, L89 

Bonnell I., Bastien P., 1992, ApJL, 401, L31 

Clarke C. J., 2009, MNRAS, 396, 1066 

Clarke C. J., Pringle J. E., 1991, MNRAS, 249, 584 

Clarke C. J., SyerD., 1996, MNRAS, 278, L23 

Forgan D., Rice K., 2009, MNRAS, in press 

Forgan D., Rice K., Stamatellos D., Whitworth A., 2009, MN- 
RAS, 394, 882 

Furlan E., McClure M., Calvet N., Hartmann L., D'Alessio P., 
Forrest W. J., Watson D. M., Uchida K. I., Sargent B., Green 
J. D., Herter T. L., 2008, ApJS, 176, 184 
Gingold R. A., Monaghan J. J., 1977, MNRAS, 181, 375 
Hartmann L., Kenyon S. J., 1996, ARA&A, 34, 207 
Herbig G. H., 1977, ApJ, 217, 693 
Herbig G. H., 2007, AJ, 133, 2679 
Ida S., Lin D. N. C, 2004, ApJ, 604, 388 
Kenyon S. J., Calvet N, Hartmann L., 1993, ApJ, 414, 676 
Kenyon S. J., Hartmann L., Hewett R., 1988, ApJ, 325, 231 
Kenyon S. J., Hartmann L. W., Strom K. M., Strom S. E., 1990, 
AJ, 99, 869 

Lin D. N. C, Faulkner J., Papaloizou J., 1985, MNRAS, 212, 105 
Lodato G., Clarke C. J., 2004, MNRAS, 353, 841 
Lodato G., Meru F, Clarke C. J., Rice W. K. M„ 2007, MNRAS, 
374, 590 

Lodato G„ Rice W. K. M., 2004, MNRAS, 351, 630 
Lucy L. B., 1977, AJ, 82, 1013 

Mayer L., Lufkin G, Quinn T., Wadsley J., 2007, ApJL, 661, L77 
Monaghan J. J., 1992, ARA&A, 30, 543 
Monaghan J. J., 2005, Rep. Prog. Phys, 68, 1703 
Pfalzner S., 2008, A&A, 492, 735 

Price D. J., 2007, Publications of the Astronomical Society of 

Australia, 24, 159 
Rice W. K. M., Armitage P. J., 2009, MNRAS, pp 709-+ 
Rodriguez L. F, Loinard L., D'Alessio P., WilnerD. J., Ho P. T. P., 

2005, ApJL, 621, L133 
Shakura N. I., Sunyaev R. A., 1973, A&A, 24, 337 
Shu F. H., Adams F. C, Lizano S., 1987, ARA&A, 25, 23 
Stamatellos D., Hubber D. A., Whitworth A. P., 2007, MNRAS, 

382, L30 

Stamatellos D., Whitworth A. P., Bisbas T., Goodwin S., 2007, 

A&A, 475, 37 
Vorobyov E. I., Basu S., 2005, ApJL, 633, L137 
Vorobyov E. I., Basu S., 2006, ApJ, 650, 956 
Vorobyov E. I., Basu S., 2008, ApJL, 676, L139 
Whitehouse S. C, Bate M. R., 2004, MNRAS, 353, 1078 
Zhu Z., Hartmann L., Gammie C, 2009, ApJ, 694, 1045 



